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Abstract 



The low energy excitation spectrum is found for a layered superconductor vortex 
with a small number of impurities inside the vortex core. All levels are found to 
be correlated. This leads to the strong enhancement of conductivity in superclean 
layered superconductors. 
PACS: 74.60.Ge,68.10.-m,05.60.+w 



1. Introduction 

In superconductors with weak pinning, the I — V characteristic displays anomalous 
propertiesB'i. Some of them are very difficult to explain in the framework of the quasi- 
classical approach. In the case of the quasiclassical approach there are three limiting cases, 
determined by values of the three parameters: the size of the gap A in the single particle 
excitation spectrum, the level spacing inside the vortex core c^o ~ A^/e^ and the electron 
mean collision time Ttr- The three limiting cases are: "dirty" limit, for tj^A ^ 1; "clean" 
limit for A ^ r^^^ ^ ujq, and superclean limit when the condition tooTtr ^ 1 is fulffiled. 

In the "dirty" limit at zero temperature the calculation of Gor'kov and Kopnini conffims 
the qualitative picture of vortex motion of Bardeen and StephenH. In accordance with the 
picture of Bardeen, Stephen the vortex core is in a "normal" state. Bardeen, Sherman! and 



Larkin, Ovchinnikovi derived conductivity in a mixed state for low temperatures and small 
magnetic field in the case of moderately clean superconductors. In this case compared to 
the previous picture a logarithmic large factor arises in conductivity. This factor is related 
to shrinkage of the vortex core at low temperatures T <^ T^. The Hall component of 
conductivity was found for moderately clean superconductors by Kopnin and Lopatini. 

The superclean case was studied in the Kopnin, Kravtsov papeiil. It was found, that 
the level spacing inside a vortex core plays the same role as the cyclotron frequency 
Uc = eH/mc in a normal metal. It was also found, that in the superclean limit the Hall 
component of the conductivity tensor is the largest one a ^ eUe/B (here He is the electron 
density in the conduction band , B is the magnitude of the magnetic field). The dissipative 
part of the conductivity tensor is smaller by the parameter (ujqt)^^. Hence the dissipative 
part of the resistance tensor is the same as in moderately clean superconductors. 

The quasiclassical approach is probably violated in the two dimensional case (in layered 
superconductors), because the excitation spectrum in the vortex core is then discrete. 

Guinea, PogorelovEl consider the dissipation in the vortex state as a result of transitions 
between unperturbed levels induced by "moving" impurities. Such a perturbation theory 
approach is valid only in the high velocity limit v ^ Vp{/S./ . 

Feigel'man, Skvorzov0 consider energy dissipation during the vortex motion as a result of 
Landau-Zener transitions between levels. They suppose, that the level distribution inside the 
vortex core obeys Wigner-Dyson statistics subject to some corrections, related to specifics of 
superconductivityEi. Such a treatment can probably be used in dirty and moderately clean 
superconductors. This method, although differing from quasiclassical approach, nevertheless 
gives for the essential range of electrical fields the same expression of conductivity as the 
quasiclassical approach. 

In this paper we consider the superclean limit. We find that in this region a new mecha- 
nism of dissipation arises. In the superclean limit no more than one impurity can be found 
at distances of order of the correlation length ^ = vp/ ^ from the vortex center . It will 
be shown that in such statistical description of level positions is impossible. If an 



impurity is placed at a distance of order of ^ from the vortex center and is weak (Born 
parameter is small), then the shift of levels is also small. It is important also, that levels 
with even and odd orbital momentum are shifted in opposite directions. The level shift 
increases as impurity comes closer to the vortex center. At some distance from the impurity 
to the vortex center levels practically cross. It is very important, that all levels with energy 
l^l <^ A cross simultaneously. 

If we neglect the weak (of order of cuq/pfC) repulsion of levels in this region, then positions 
of levels as a function of the distance from the vortex core to the impurity form two families 
of crossing straight lines. Outside of the dangerous level crossing region these lines are 
practically horizontal (see fig.l). The size of the dangerous zone, where the level hues can 
be considered as crossing, depends on the Landau-Zener parameter and hence on the vortex 
velocity. For the vortex velocity V in the range ujq ^ ppV ^ LUolA/ep) the distribution 
function of excitations inside the core of the vortex does not change until the impurity 
comes into the dangerous zone. But it changes essentially when impurity goes through the 
dangerous zone. Excitations, that arise when impurity goes through this zone, determine 
the value of the dissipative part of conductivity. Such mechanism of dissipation is essential 
for the magnitude of the electrical field E lying in the range 



Here c is the velocity of light. As we prove below, this mechanism of dissipation leads to 
the dissipative part of conductivity axx being equal to 



where Qq is the distance from the "dangerous" region to the vortex center oq ~ ^ is the 
Born parameter which is equal to the phase shift of an electron scattering off the impurity. 
Usually its value is of order of one ^ ~ 1. Hence parameter (pf^o) is much larger than 
one and axx essentially exceeds the value, obtained in the framework of the quasiclassical 
approximation. In the range ppo, ^ ooQTtr ^ 1 the Hall angle is small. Nevertheless the 
dependence of (Txx on the scattering time r^j. is the same as at large values of the Hall angle. 




3 



2. The low energy excitations spectrum for an impurity at the distance a from the 
vortex center in the range a » ^(A/ei?)^/^ 



The excitations spectrum E in the vortex state can be found as a solution of the eigenvalue 
problem for the system of equationslllll 



yf2j 



(2) 



-2Wia^-f^ + yi^)-E; A(r) 

where A is the order parameter, fi is the chemical potential (Fermi energy fi = ep), V{r) is 
the potential of impurities. We suppose here, that the magnetic field B is weak {B ^ Hc2) 
and omit the vector potential in the eq.(2). Below we consider the two-dimensional case. 
We suppose also, that there is only one short range impurity (with the interaction radius of 
order of pp^) inside the vortex core. 

In our problem the order parameter A in the absence of the impurity is given by the 
expression 

A(r) = A(r)exp(iiv?) (3) 

where (f is polar angle, r = |r|. The low energy excitation spectrum in the absence of the 
impurity was found in the papeilll. The system (2) possesses a very important property: if 
E is an eigenvalue with the eigenfunction (/i,/2), then —E is also an eigenvalue and the 
corresponding eigenfunction is (/g , — /i )• This property holds in a magnetic field too. 
The low energy excitation spectrum E^ is given by the equationlll 

E° = -in - 1/2)^0 (4) 
where 

JO ppr Jo 

K{r) = r driA{ri)/vF ■ (5) 
Jo 



If Kramer-Pesh effect takes place, then with a logarithmic accuracy we obtain from the 
equations (4), (5) 



/A 

c.o = -ln(- 



A = A 



(oo) 



(6) 



The eigenf unction, corresponding to the eigenvalue (4) is 





[/.] 


/ 


fn — 








V 



(7) 



where c is the normalization constant, J„(x) is a Bessel function, n — 0, ±1, ±2, .... 

For the excitation spectrum E in the presence of an impurity inside the vortex core, we 
obtain from the eq.(2) the following system of equations 



det ((£ - £;) + i) = 



(8) 



where the operator A is given by its matrix elements. In the basis (7) we have 
^ V(r-ai): ^ _ 

fn> , (9) 



^kn —< fk 



V"(r-a); 
V 0; -^(r - a) ) 

a is the position of the impurity relative to the vortex center, and 

r — A 

tfcn "kn-'-'n • 



(10) 



In the eq.(9) essential are values of r such that r 3> Pp^. So we can use an asymptotic 
expansion of Bessel functions to find matrix elements A^n- A simple calculation gives 



A,^ = e^^^-")^" |/i(a) cos (^^^t±3)j _ /^(«) 



'n{n + k) 



(11) 



where is the polar angle of the vector a, and the quantities /i^2 are given by the equation 

cos(2p^(a+^)) ) 



^sin(2p^(a + ^)) ^ 



(12) 



In the eq.(12) the normalization constant C is equal to 
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= {2n /"dre-2^('^)}-i (13) 

If there are several impurities inside the vortex core, then the operator A in the eq.(8) is a 
simple sum over all impurities. Therefore 

At^-T^AUai) (14) 

i 

where Akn{ai) is given by the eq.(ll). It follows from the equations (11), (14), that the 
transition matrix elements A^n are separable. That is A^n can be presented as a finite sum 
of terms of the type A^Bl 

j 

As a result we can obtain an expression for the excitation spectrum in an explicit form. If 
only one impurity is placed inside the vortex core, then we obtain from the equations (8), 
(11), (15) the following equation for the excitation spectrum 



det 



/ N+l ^ N+l \ 

1 + -^1 E cot+e'^ -^2 E . B 

L=-N L=-N '^^^^^ 

N+l ^ N+l 

h E , .r; 1 - -^1 E 



(16) 



For the linear spectrum given by the eq.(4), we obtain in the limit N — > +oo 

^1 1 TT / /I E 
> = Ctg TT - H 

£ ^ "i^'^ (i " ^)) ^''^ 

With a help of the eq.(17) we reduce the eq.(16) to the form 

l + ^l(M^ + ^^^0 (18) 

4^0 ^ocos(5) 

It follows from the eq.(18), that the low energy excitation spectrum is strong correlated even 
in the presence of an impurity inside the vortex core. If Eq is a spectrum point, that is if 
Eq is some solution of the equation for the spectrum (18), then all solutions of the eq.(18) 
are given by the equation 

6 



E=±Eo + 2uJoN, N = 0,±l,±2. (19) 

Hence the discrete spectrum is given by two sets of equidistant points. 

Functions Ji^2 are periodic with the period tt/pf- Both are defined by the same function 
with shift by a quarter of the period. The amphtude Ii,am of these functions is a smooth 
function of the parameter {a/^) and is given by the eq.(12). 

With the accuracy of ujq^A/ef) a point exists such that 

h.aM = -— ; hiao) = . (20) 

71 

Hence at the points + 6a given by equation 

Sa = (^-^^ N , = 0, ±1, ±2... (21) 
we have 

- T (^) (-) 

Equation (22) means, that in a vicinity of the points of the trajectory of the vortex, given 
by the eq. (20), there is a set of points, separated by the distance 6a, where spectrum fines 
are practicaUy crossing (see Fig.l. In Fig.l tfie quantity 6a = a — a^ is sfiown, wfiere ao 
is given by tfie eq.(PD|).) Tfie vicinity of sucfi points we denote as lagoon. If tfie impact 



parameter of tfie trajectory is smaller tfian some critical value, tfien on sucfi a trajectory 
tfiere are two lagoons. Wfien tfie vortex moves tfirougfi tfiese two lagoons many excitations 
are created inside tfie vortex core. Tfie contribution of tfiese excitations to tfie dissipative 
part of conductivity will be found below. 

Tfie situation of fiaving several impurities inside tfie vortex core is considered in Ap- 
pendix. We prove tfiere, tfiat tfie strong correlation in tfie level positions survives for two 
impurities inside tfie vortex core. We can make a conjecture tfiat tfie strong correlation in 
level positions exists for many impurities inside tfie vortex core, wfiile tfie condition It^ ^ C 
if fulfilled. But on a vortex trajectory of a general status lagoons witfi large number of 
practically crossing energy level lines do not exist, if tfiere are two or more impurities inside 
tfie vortex core. 



3. One impurity at small distances (a <^ ^(A/e^?)^/^) from the vortex center 



First of all we shall consider one impurity with a short range (of order of pp^) potential 
placed exactly at the vortex center. 

At distances p ^ pp^ from the vortex center we can use for the solution of the eq. (2) the 
quasiclassical approximation with the first order correction terms. Indeed these correction 
terms will give an expression for spectrum. And at small distances of order of pp^ we can 
omit nondiagonal elements in the eq.(2). As a result we obtain the following expression for 
the spectrum 

On-l — On 



-{n - l/2)ujQ + vpt^ 



POO 

I \ dpe~ 
Jo 



■2K{p) 



(23) 



where On is the scattering phase in a state with angular momentum n in the presence of the 
impurity potential ^(r). The corresponding eigenfunction is given by the expression 



\f2j 



c 



-K{p) 



e*"'^( J„(PfP + On) + Jn{PFp + On-l) 
-e*("-l)'^(J„_l(p^p + On) + Jn-1{PFP + ^n-l)) 



(24) 



where C is a normalization constant. Suppose now, that the impurity is placed on a dis- 
tance a from vortex center, such that |a| ^ ^(A/e^)^/^. In the equation (2) we make a 
transformation to the coordinate system with the origin at the impurity. Then we obtain 



/ 



A^I^-p + V{v)-E- 



|A|e*^+ (af ) (|A|e^'^) \ ( f:^ 



|A|e-^+(a|:)(|A|e- 



-i<p\. 1 (9^ 
> ' 2m ar2 



+ p-V{v) 



E 



\f2j 







(25) 



From the eq.(E5|) we obtain the following equation for the excitation spectrum 



det Ui - E)+A 







(26) 



where e^n = EnSnk and En is given by the eq. 
elements A^n in the basis defined in the eq.(E3). 



A 



kn 



--< fk 



0; 



(a|:)(|A|e^^; 



a#)(|A|e-^^); 



fn > 



The operator A is given by matrix 



(27) 



A simple straightforward calculation making use of the equations 



gives 



Jo p I 



5k,n-ie''^^ COS } 



e cos + 



(28) 



In the eq.(28) the constant C is given by the eq.(13). It follows from the eq . (pSQ , that in 
the operator in the eq . (p6D , only diagonal and near diagonal elements are nonzero. Now we 
define the function B[I, E, n) in the following manner 

J2 pQg2 ( ^n-l-^n + l ^ 



B{I,E,n-l) = -En-E 
where 



B{I,E,n) 



(29) 



roo A 

I = naC^ / dp^-^ ■ e-^^^p\ I/uq = pFa/2. 
Jo p 



With the help of the function B we reduce the eq.(P^) for the spectrum to the following 
simple form 



det 



B{I, E, 1); -Je*^" cos {^^) ; 0; 



60-02 



0; -Je-^^^ El - E; - Je^^" cos (- 
0; 0; Je-^"^" cos (^a^) ; 5(/, -E, 1) j 

The equation (|29|) means, that B{I, E, 1) can be presented as an infinite fraction 

/2cOs2((^^l-03)/2) 



(30) 



B{I,E,1) = -E2-E- 



-E. - E - 



/2cos2((^^2-^4)/2) 



(31) 



-E.-E- 



Pcos\{e,-e,)/2) 



. B{I,E,n + l) 

The fraction (^1]) converges very quickly, if for B{I, E,n + 1) we use the expression 



B{I,E,n+l) 



n + l/2)o^o -E + d{{n + l/2)u;o - Ef - AP 



/2 



(32) 



{n + l/2)cuo ± ^ > |/ 
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(33) 



Suppose now, that the impurity is of a small size, so only 5'-scattering is essential. Suppose 
also, that the impurity potential is of order of the atomic one and hence the inequality takes 
place 

eo = Atg(0o/2) > ivo (34) 

Then in the first approximation the equation (pOD for low energy excitations uncouples to 
two independent branches 

5(7,^,1) = and B{I,-E,1) = (35) 

And hence we obtain two independent families of spectrum lines. Of course in this 
approximation they will cross. And only in the next approximation with respect to the 
parameter (uo/eo)^ a gap in the crossing points will open. 

For small values of E we have 

B{I,E,1) = B + aE (36) 

where 

dB{I,E,l) 
" = dE 



Inserting the expression (^) into the eq.(|5DD we obtain the following equation for the lowest 
energy level near the crossing points 

E^a^l = B\lcos{eo/2)f + {e^B + (Jcos(0o/2))')' (38) 
Hence the value of the gap 5 near this crossing point is equal to 

i = l^"-'"'/^'!' (39) 

\ae„\Jel + {Imi{eol2)f 



In the eq. (pSD the quantity a should be taken at the point 
dB{I,E,l 



a 



dE 



2^ EO(Jcos(eo/2))2 (40) 



(i-cos(eo/2)) 



2T7I 
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By the order of magnitude the value of the gap S is given by the equation 

5 ~ A(a/0' (41) 

It is possible to keep in the expansion of the order parameter A with respect to the 
shift a in the eq.(^) terms up to the second order in a. Then the operator A will have the 
following nonzero matrix elements: Ann] ^n,n±i; ^n,n±2- The equation p6|) in this approach 
enables us to determine the excitation spectrum up to the shift a of order of 

a ~ aNepy/' (42) 

At the boundary of this region the gap S, given by the eq . (|4l|) , is of order of ujq. Hence 
energy levels can "cross" only in the lagoons (eq.(22)), or when the impurity is placed near 
the vortex core (eq.(p^)). 

The equation (pOD for the spectrum can be reduced to the form 



(^11 cos' |^ ^o-^2 ^ + BiI,E,l) + BiI,-E,l) J 
' B{I,E,1)-B{I,-E,1) _EP^^^, I^ 9o-92 ^J ^ 



-B{I, E, l)BiI, -E, 1) + cos^ =0 (43) 



Ei ^ ' ' ^ ^ ' ' ^ ' V ^0 V 2 _ 

Near the crossing points the last two terms are small (of order of {I^ / EqY ) and lead to 
the repulsion of spectrum lines. In zero approximation neglecting two last terms the eq.(^3|) 
uncouples to give two families of independent spectrum lines. 

Energy levels as a function of the shift a (or quantity / ~ a) are given in the Fig. 2. 
Inside the circle in the Fig. 2 the equations for two spectrum branches are 

E/ojQ - 1.073 = 5(0.1335t ± Vl + t^); (44) 
t = 1.343(//u;o - 2.3172)/(5; (45) 
6 = 3.2- 10-^ (46) 

The contribution of these "crossing" points to the dissipative part of conductivity will 
be discussed below. 
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4. Landau-Zener tunnelling near the crossing point of spectrum lines 



Near the crossing point of spectrum lines, given by the eq.(18), we can put 

^ = -l + y- ^ = 2pfX- X = Vt (47) 

Where V is the velocity of the vortex. For two close spectrum points E± we obtain from 
the eq.(18) the following value 



E± = ±£(,)^ , = ^y^ + {2pFXf (48) 

The usual Landau-Zener consideration leads to the following value for the probability 
for a "particle" to remain on the same branch after collision 

With a help of the eq.(^) we will find the energy, transmitted to the vortex at one collision 
with impurity. 

5. Dissipative part of conductivity in superconductors 

Let us first consider the superclean limit, when only one impurity can be found inside the 
vortex core. Then for small values of the vortex velocity V , such that 

V < VFi^/SF? , (50) 

an excitation can arise only at spectrum line crossing points. Such crossing points are located 
only in lagoons (eq.(20), (21)), or if the vortex center is close to the impurity (eq. (p9|)). 

Suppose, that between two consequent collisions the system goes to the equilibrium due 
to inelastic scattering processes. In such a case the dissipative part of conductivity is directly 
connected to the energy, stored by vortex after one collision with impurity. 

Consider first processes happening in lagoons. Suppose, that po is the impact parameter. 
Then near the point (eq. (20)) we have 
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If-'-^^' ^-i(^) .VT^W^ (51) 



ao 



Periodicity 5x of a crossing point with respect to x is 

5x = (52) 

Pf^/I - (po/flo)^ 

Hence the full number of crossing points in a lagoon, with effective transition of "particle" 
to the other branch, is 

~S^~ ~^ TK^) ^^^^ 

V dp 

To obtain the eq.(^) we use the eq. f^Ql) . We get with a help of the eq.(p3D, that the vortex, 
after passing through two lagoons (one collision with impurity), stores the energy 6E being 
equal to 

{dIi,am/dp)lo 

As a result, we can estimate the contribution of transitions in lagoons to the dissipative part 
of conductivity cr^]^ 

(1) _ 256 apnimp (^jpp . . 



.am 



ao 



where (j)Q = vr/e - flux quantum, B - magnetic field value. Consider now the energy, dissipated 
as an impurity passes near the vortex core. From the eq.(0), (|39D we obtain the following 
expression for the excitation spectrum near the crossing points 

/ j6 \l/2 

E = ±{-^ + {uoPF{5a)/2f (56) 



where 



dB{I,E,l) 

72 = -Q-j ; I = uJoPFa/2 (57) 

{6a) is a shift from a "crossing point", the parameter 72 is of order of unity. According 
to the eq.(^), the probability for the particle to remain on the same branch after a 
collision is equal to 
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Hence the number of excitations N, that arise in the vortex core, when the impurity 
passes through the vortex near its center, is 

""--r^ U(-o/>./2)^j ^''^ 

These transitions give the following contribution to the energy dissipation 6E per volume 
and time unit 

SE = 2aVuoN^ni^p{B/<Po) (60) 
Equations (^5]), (|60|) completely determine the dissipative part of the current 

^''---^>.^''+^- — [,(^„p,/2).) (^/B) (61) 
By the order of magnitude we have 
/ dli.am \ 2a;o 



dp I vrf 

/ ao ^ 



(62) 



Hence, in the range of velocities V, such that 

V/VF > [A/bpY (63) 

the second term in the eg .(1611) is smaller than the first one. By the large parameter 

^-^^{ep/Af (64) 

the dissipative part of conductivity, given in the eq.(0), exceeds the quasiclassical value for 
conductivity in the two dimensional case. The expression (p5D for conductivity was obtained 
under the assumption, that the number of crossing point in a lagoon eq.(|53D is large. This 
condition gives the same restriction on the value of the velocity as given in the eq.(B3). 
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6. Conclusion 



In the two dimensional case the excitation spectrum in vortex core is discrete. This results 
the strong increase of the dissipative part of conductivity Or^^ (eq-(l)) comparing to its 
value obtained by the quasiclassical method. It is very probable, that the strong increase of 
conductivity On^x at low temperatures, obtained in experimental paperJ^'i, is related to the 
phenomena, considered in this paper. For the detailed comparison to the experiment both 
experimental and theoretical investigations are necessary. 

Consider now the applicability region of the equations. (^^D, (pT]), (|5^). We do not see 
any restriction for the temperature T in the framework of the model used except T -C T^. 
In the range T ^ excitations with the energy e ~ A are essential. In this region the 
excitation spectrum is not equidistant. Probably its more essential, that in our model the 
main mechanism of dissipation is the scattering of excitations off impurities. In real high-T^ 
compounds such an assumption holds only in the low temperature region. The predicted 
effect should vanish in the strong magnetic field region H ~ iJc2- 

More complicated are restrictions for the electrical field magnitude or for the vortex 
velocity. The restriction vpi^l^j Spf" <^V <^ vpi^t^j spf' ■, given in the introduction, holds in 
our approximation. For the vortex velocity V ^ vp{^/£pY the adiabatic consideration is 
inapplicable. An impurity on the distances a < ^ from vortex center leads to the transitions 
of quasiparticles to highly exited states. In the region V '^VFi^/EpY ^^e small gap in the 
excitation spectrum of order of ujQ{/S./eFY leads to the adiabaticity of motion in lagoons. In 
this case the impurities, that pass at small distances of order of ^(A/e^?)^/^ from the vortex 
center (second term in the eq.( ^l])) give the main contribution to the energy dissipation. It 
is essential, that even in the frame of our approach the I — V characteristic changes its form 
for different ranges of the velocity. 

There are also some other mechanisms, that can change the I — V characteristic. At 
small values of vortex velocities pinning can start to be essential. It can reduce dissipation. 
The pinning force is strongly dependent on the interaction between vortices and hence on 
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the magnetic field magnitude and anisotropy factor. In superclean superconductors pinning 
can be considered weak. 

Tfie tunnelling of excitations to the neighbouring planes can reduce dissipation at large 
velocity. If the anisotropy factor e is small then the spectrum near the crossing points does 
not change much and the Landau- Zener tunnelling probability does not change. But near 
the points, where the energy of a quasiparticle is close to the unperturbed position of a level 
from the neighbouring plane, a notable probability arises for tunnelling of a quasiparticle to 
the neighbouring plane. As a result, the quasiparticle energy does not increase any more and 
dissipation will be smaller than given by the eg . (|5^ . The energy relaxation of quasiparticles 
can change I — V characteristic as for large values of vortex velocity so at small one. If the 
energy relaxation time is large enough, then effective heating of excitations inside of the 
vortex core takes place. The effective temperature reaches a value of order of The 
value of axx decreases in this case. On the other hand, strong energy relaxation can lead 
to the equilibrium inside of one act of tunnelling. This phenomenon will decrease a^x value 
in the range of small vortex velocities. Now the physical reason of the energy relaxation 
in high-Tc superconductors at low temperatures is not clear. And we can not make any 
quantitative estimation of this effect. Equations (|5BD, (^Tj) are valid in the superclean limit, 
if at distances of order of ^ from the vortex center there is no more than one impurity. If at 
a distance of order of ^ from the vortex center there are two impurities, then the probability 
of a level crossing decreases strongly. As a result, the value of axx decreases too. 

(7 XX ^ 

function of impurity concentration has a maximum at ujQTtr of order of one. 
The question about the value of the impurity concentration, for which the matrix can be 
considered as random demands an additional investigation. 
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Appendix A: A large number of impurities inside of the vortex core 



If inside the vortex core a large number of impurities is placed, then transition matrix 
elements A^ni^i} {^i is the position of i-th impurity), are given by the eq.(14). Equation 
for the spectrum in such a case is 

det(i + C) = (Al) 
with the matrix elements of the operator C being equal to 



If M impurities are placed inside the vortex core, then the size of the matrix in the eq. (|Al|) 
is (4Mx4M). Due to symmetry properties of the elements Aj^n (eQ-(ll))) this matrix can 
be easily reduced to the size (2Mx2M). The structure of this matrix M is simple. On the 
diagonal are placed blocks (2x2), defined only by one impurity. The second type of blocks 
are blocks (2x2), that give the interference contribution of a pair of impurities (ij). 

To clarify the structure of the matrix M in the general case, we give below the explicit 
expression of the matrix M for two impurities inside the vortex core. 



M 



(A3) 



1 + S^tg (f + ' S^tg (f + if) ' 2!^! ' 2!^! 

Jj_^*. J2_^*. 1 _|_ ZLSpfrr r 2L _|_ 2LR\ ■ ^rfp- (- 4- —\ 
2c^o 1' 2<^o 1' ^ 2<^o g U ^ 2uJoJ ' 2<^o g U ^ 2lvqJ 

\ 2uJo ' 2^0 ' 2(^0 2(^0/ ' 2t<;(,^^6 V4 2uJo J / 

In the eq. ( |A3| ) the upper index in Ij means the number of impurities (j = 1, ...M) and the 
lower index i can be 1 or 2. The quantities z, Zi are defined by the equations 

^2iL{ipa^-Lpa2) 



E 

L=—oo 



L-1/A + E/2uoq 



00 gi(2L + l)((/Pai-'/'a2) 



From the eq . ( |A4|) it follows, that 
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zi{E) = -^*(-E)e^('^-i-'^»2) (A5) 
A straightforward calculation gives for the quantity z{E) the following expression 



z{E)- 

where 
Y 



TT 



sm 



7r TvE \ 

,4 2ujoJ 



-i7r(l/4-£;/2wo)+jXaia2(l/2--B/aJo) 



(A6) 



iVai - </?a2)/mod TT > 



(A7) 



From the eg . ( |A3|) we obtain 

nil 



detM 



1 + Ji + 



Uq cos I — 



1 + J2 + 



Uq cos 



(-) 



+ 



T'^T^ I rr M 

-(ziz + + -^^r^ Ji J2- 



L 



l£i|_ 



ji + ^j,ctg - + — 

TT V 4 2uJq 



. 2ulIq ^ (tx txE 
jJ + ^Jictg -- — 

TT V 4 2(jJo 



2LUn ^ / TT TvE 

— ^2Ctg - + — 

TT V 4 2a;o 



2uJn 



TT 







(l 


2cuo/J 



(A8) 



where 



(A9) 







Note, that the quantity zzl is pure imaginary. Hence the only term in the eq. that 
has a periodical dependence with respect to the energy E with the period, differing from 
2ljq, drops out. As a result, the expression for the spectrum (19) holds also in the case, if 
there are two impurities inside the vortex. Our conjecture is: for the low energy excitation 
spectrum eq.(19) is correct in the clean limit (rA ^ 1), even if there are a lot of impurities 
inside the vortex core. 

With a help of the equations (|A5|), (|AB|), ( |A^ ) we obtain the following equation for the 
excitation spectrum 




detM = 0, 



det M = 1 + Ji 



TT^ ( . 2UJQ ^ fTT 7r£'\\ / „ 2uJo ^ fTT TlE 



For the I — V characteristic in the mixed state decisive is the value (or the existence) of 
the gap in the excitation spectrum. For this reason we will calculate the value of det M for 
E = {). Simple calculations with a help of the eq.( [Al(j| ) give 

det M(E = 0) = ^(1 - Ji) + ^(1 - h) + (1 + Ji)(l + J2) (All) 
The expression ( |Allj ) is nonnegative. It can be written as a sum of nonnegative terms 



det M{E = 0)= 1 + ^ + ^ 

\ 2ujq 2ujq 2ujq 2ujq 



2 

I + 



It means, that on the vortex trajectory of a general status, lagoons with large number 
of practically crossing energy levels lines do not exist. Now we can formulate the following 
conjecture. In the clean limit (rA ^ 1), if several impurities are placed inside the vortex 
core , the quantity det M{E = 0) is a sum of nonnegative terms of the form ( |A12|) . 

Now, for small values of <^ ujq, we obtain from the eq.( |Al(]| ) 



detM(E) = detM(E = 0) 



^2J,J, - ^(1 - J,) - gL(l - Ji)| (A13) 
The equations (|A11|) , (|A12|) , (|A13D completely determine the excitation spectrum near the 



crossing points of the energy levels. 
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FIGURES 

Fig. 1. The excitation spectrum as a function of the impurity distance from the vortex center. 

The parameter {-k^ /2L0Qpp){dIi^am/d(i ■ ■ •) equals 0.02. The quantity 5a = a — oq, with oq is given 
bytheeq. (20). 

Fig. 2. The low energy excitation spectrum at small distances a from the impurity to the vortex 
center; I/uJo =pFa/2, see eq.(29); 0^1; Eq/uq = 50. 
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